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TRANSPORT MODEL OF NUCLEON-NUCLEUS REACTION 
Abstract 

A simplified model of nucleon-nucleus reaction is developed and some 
of its properties are examined. Comparisons with proton production 
measured for targets of 27 hi> 58 Ni , 90 Zr, and 209 Bi show some hope for 
developing an accurate model for these complex reactions. It is suggested 
that binding effects are the next step required for further development. 
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The understanding of the interaction of energetic charged particles 
with matter is of importance to radiation protection in space as well as in 
astrophysics and radiotherapy (ref. 1). Fundamental to proton transport is 
the understanding of the nuclear reactions involved. Considered herein is 
a simplified model of the nuclear induced reactions with the intent of 
developing a simple formalism which is accurate enough for transport 
calculations. 

Theory 

Many years ago, Serber proposed that high-energy nuclear reactions 
consist of two distinct phases when viewed on a time scale (ref. 2) First, 
the passage of the initiating event is on the order of the transit time 
across the nuclear volume requiring only 

tj « 5 fm/c = 1.7 x IQ" 23 sec 

and is short compared to the internal motion within the nuclear interior 
(ref. 2). The second phase of the reaction, which involves the residual 
excitation of the target nucleus after the passage of the projectile, is 
characteristic of the internal nuclear structure and occurs on a much 
longer time scale (~10“ 15 sec). 

To treat the first phase, an approximate Boltzmann equation will be 
applied to nuclear matter (refs. 3 and 4). Particle spectra at the end of 
the first step will be evaluated. These may be used in a compound nucleus 
model to further define the Serber second phase. 

The transport equation in nuclear matter is 

Fn»v + P CTlip(x,n,E) = fpcrf (E ,n ,E 1 ,n 1 ) 4 » (x,fi 1 ,E ' )dn 'dE ' (1) 
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where a is the two-body cross section for interaction of the field i{> with 
individual constituents of the nuclear matter, p is the constituent 
density, f(E,n,E',n') is the secondary production spectra in the 
interaction events, and iji(x,fi,E) is the field function of the transported 
particles. 

Past studies have shown that the straight ahead approximation is fully 
justified in space radiation shielding applications (ref. 1). This 
encourages us to ignore the angular dependence in the nuclear problem as it 
is ignored in the transport problem. All particles are then assumed to be 
produced in the forward direction as would be expected at high energies due 
to relativistic effects as well as the dynamical nature of the two-body 
force. This causes simplifications in equation (1) as 

rfz + pol«fp(b,Z,E) = /°°pf(E,E l )ij)(b,Z,E' )dE' (2) 

A 

where the direction of motion is along +Z, and b is the impact parameter 
vector. The integral form of equation (2) may be written as 

Z + 

* -f p(b,y)cr dy + 

«K(b,Z,E) = e ^(b.E) 

V + 

z f P(b,y)dy + 

+ f e p(b,Z' )dZ' / f (E ,E 1 )i(i ( b,Z,E 1 )dE‘ (3) 

-oo E 

where ip -f ( b ,E ) is the incident particle fluence at impact parameter b, and 
its coefficient in equation (3) is the appropriate attenuation factor. In 
what follows, the incident fluence is taken as 

* i (b,E) = «(E-E 0 ) 

corresponding to a unit fluence with energy E 0 . 


( 4 ) 



Attenuation of the Primary Beam 
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The attenuation of the primary beam is given by the inhomogenous term 
in equation (3) as 

Z 

„ -J o P (b,y)dy + 

♦ 0 (b,Z,E) = e ♦ i (b,E) (5) 

which are those particles passing through to a position Z without 
collision. The number which suffer collision before reaching Z are 

Z -*■ 

-f crp(b,y)dy + 

) ^ i (b s E) (6) 

from which the absorption cross section is obtained by integrating the 
absorption factor over all impact parameters as 


abs 


= 2ir f (1 - e 


oo ->■ 

\f crp(b,y)dy 

— oo 


)bdb 


(7) 


in the limit as Z becomes large. Evaluating for a uniform nuclear density 
yields 


abs 


= 2ir f (1 - e 0 


-2pa/R 2 - b 2 


)bdb 


= ttR 


0 

2 2tt 


(2p 0 a) 


2 


(8) 


where p 0 is the nucleon density, R is the nuclear radius, and a n (Z) is 
related to the exponential integrals. In the limit as p 0 a->-«> yields 
,2 


“abs * ” R 


(9) 


as represented by an absorbing disk, and in the limit as p 0 ct>o yields 


a abs “ " t Po 0R ) 


( 10 ) 
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where the first term results from scattering from each individual 
constituent and the second term is the first shadow correction. 


Two-body Scattering Cross Sections 


The nucleon-nucleon differential cross section is represented as 

W ■ Ae +St 


(ID 


where b is the usual slope parameter and t is the square of the 
four-momentum transfer. A is related to the total cross section and to a, 
the ratio of the real -to-imagi nary part of the scattering amplitude 
(ref. 5) as 


A = (1 + « 2 ) (12) 

The factor of 2 in equation (12) approximately accounts for the backward 
peak in the cross section. The energy transfer cross section is then 


# = 2mc 2 Ae- 2mc2f?(E o " E) 


since 


(13) 


t = -2mc 2 (E 0 - E) (14) 

where m is the rest mass of a nucleon and (E 0 - E) the energy transfer. 
Considering the target nucleon at rest yields the struck nucleon recoil 
spectrum as 



= 2mc 2 A 


-2mc 2 BE 


(15) 


Thus the spectrum of scattered particles is given by 


f(E.E') = 2mc2Are-2mc 2 B(E' - E) + e -2mc 2 BE] 


(16) 
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for use in equation (3). Since the struck nucleon is immersed in a "sea" 
of nucleons, not all energy transfers are allowed. Only those for which 
the struck particle is lifted above the fermi surface may scatter. The 
appropriate cross section is then 

9 E o -2mc 2 3(E n - E) 
o(E 0 ) = 2mc A f e dE 

= f a 0 f2mc 2 8(E 0 - E p )) (17) 

where the contribution from Ep = o is the Pauli blocking factor. 

Multiple Scattering Series 

Now consider the transport equation neglecting the recoiling target 
particles 


f| 7 + p(b,Z)al*(b,Z,E) = p{b,Z) / 


2mc 2 A 


; -2mc 2 3(E' - E )^( baZaE > ) dE * 


(18) 


This is solved by treating the left side of equation (18) as a perturbation 
with the first term 

Z -v 

+ - f crp(b,y)dy 

* 0 (b,Z.E) = e ■" ♦ 1 (b.E) (19) 

and assuming boundary condition in equation (4). The first perturbation 
term is 


- f a p(b,y)dy 


[fy + p(b,Z)a] ^^(b.Z.E) = p(b,Z) e 


2mc Ae 


-2mcS(E n - E) 


(20) 
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for which 


z -f «p(bjf)<|y „ -2mc 2 S (E -E) 

^(b.Z.E) = / e p(b,Z' )2mc 2 A e 0 dz’ 


( 21 ) 


The second perturbation is then 

Z + 

+ l -/ op(b,y)dy + 

4» 2 (M,E) = f e P (b,Z' ) 


v + o o -2mc Z 8(E n - E) 

x f p(b,Z' )dZ‘ * dZ 1 (2mc A) 2 (E q - E) e 0 


( 22 ) 


and the third iteration as 


* 3 (b,Z,E) 


Z + 

l ‘f "Ofbjr)* „ z . 2 „ 

f e "" p(b,z')f p(b,Z")f c(b,Z , ")dZ , dZ"dZ"' 


1 2 3 2 _ 2mc 8(E - E) 

x j (2mc 2 Ar (E 0 - E) 2 e 0 


(23) 


The general term may be written as 


* n (b,Z,E) = <j> n (b,Z) F n (E) 


(24) 


where 


P n (Z + /R 2 - b 2 ) m -2p (Z + /R 2 + b 2 ) 

2 rt ^ 


2 , u2> 


<t> n ( b > z ) = — r(n + 1) e 


(25) 
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The F n (E) is 


2 n < E o " E > 
F n (E) = (2mc 2 A)" 


n-1 


-2mc^(E 0 - E) 


( 26 ) 


The even smaller perturbations from the target constituent recoil spectrum 
will now be treated. 


Multiple Production Series 

The target recoil constituents are of inferior energy compared to the 
incident beam and are assumed to be reabsorbed by the surrounding nuclear 
material if they further interact after initial production. For each of 
the above terms there is a production term, the first of which is given by 

Z -v 

+ + + -f crp(b,y)dy 2 

^ + p(b,Z)al *{(b,Z,E) = P (b,Z), e 2mc 2 Ae~ 2mc pE (27) 

The solution is 


*J(b,Z,E) = ^(b.Z) 2mc 2 A e' 2 " 10 ^ 


(28) 


= 4> 1 (b,Z) F 1 '(E) 


Similarly 


*;(b,Z,E) 


* n (b,Z) f;(e) 


(29) 


where 


r , / r \ _ -2mc 2 6E r A^ n a n-2^ 2mc b(E o " 

F n (E) 2mc A e ( ? ) 


( 30 ) 



Multiparticle Production 
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All particles escaping the nuclear reaction may now be 
summing terms and allowing Z to approach infinity as 

^(b,E) = ^(b.E) + tf» 2 (b,E) + ^(b.E) + . . . 

■> + 4 * 

+ 1 (b,E) + tp 2 '(b,E) + if> 3 ' (b,E) + . . . 

The particle production cross section is then given by 


da 


= 2tt f i|»(b,E)b db 


dE 


= y or i rF-(E) + F 1 , (E)1 


where the for a uniform nuclear model is 


* R -*■ 

<x.j = 2 it f <j> - ( b ) b db 
o 


2, P Q 1 2 1 


r(i + l)(2p Q a) 


Z TT a i + i ( 2 p 0 ctR ) 


The total production cross section is found by integrating 
to obtain 


A A 


calculated by 

(31) 

(32) 

(33) 

(32) over energy 


Vod = \ CT i rF i + F i ' 1 


(34) 
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where 


F. 


i ■ r F - (E )dE 
E F 


. ,V °1 - l f2mc e ( E Q- E F» 
6 J r(i ) 


and 


F T 


/ F. * {E)dE 
E F 


= F, e" 2,nc bE F 


It follows from conservation of particles that 


( 35 ) 


( 36 ) 


A A 

T < ’i F i = °abs 

The multiplicity of the transport event as herein approximated is 
m - CT prod^ a abs 


-2mc^BE F 

= 2 + (1 - e h ) (37) 

which is an indication that the recoil particles need to be treated in more 
detail since intranuclear cascade codes indicate 2 m < 2.5. 



Proton Production in Nuclear Reaction 
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It may be easily shown that the proton production spectrum is 
approximated by 

»it F i< E > ♦¥ E >i < 38 > 

where Z is the atomic number and op the appropriate proton cross section 
for the projectile. Results of equation (38) are compared with the 
experiments of W u et al. in figure 1. Quite good agreement is obtained 
at the highest energies (>40 MeV). The low energy portion of the spectrum 
is adversely affected by high order multiple scattering terms adding to the 
production of protons. These terms would be limited in their contributions 
if the bulk binding potential are properly treated. A more complete theory 
should improve the agreement. 
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An Appendix of Useful Functions 
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Several functions that are closely related to exponential integrals 
and gamma functions are herein defined. 

- Z 

o (Z) = f t n e" r dt (Al) 

o 

with special value 

o 0 (Z) = 1 - e“ Z (A2) 

and recurrence relation 

“n (Z) = n “n-l (Z) “ Z " e ‘ Z (A3) 

A useful integral is 

yz) = ; o„(t)dt (A4) 

0 

with special value 

A q (Z) = Z * 0 (Z) (A5) 

and recurrence relation 

A n (Z) = nA n - 1 (Z) " a n ( z ) ( A6 ) 
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Useful functions related to a n (Z) are 

B n (Z) = / Z e t a n (t)dt 

o 

1 r 

c (Z) = r tea (t)dt 
o 

with special values 

B q (Z) = e Z - (1+Z) 


C 0 (Z) = (Z-l) e Z + l - Z 2 /2 
and recurrence relations 

B n (Z) = nB n _ X (Z) - Z" + V(n + 1) 

C n (Z) = nC n _ j(Z) - Z n + 2 /(n + 2) 


(A7) 

(A8) 

(A9) 

(A10) 

(All) 

(A12) 
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ENERGY (MeV) 


Figure 1-Protm spectra from proton impact on nuclei at 90 MeV. 
Short dash is the calculations of Wu et al . . long dash is the 
experiments of Wu et al. and the full line is the transport 
theory. 
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